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The QCD sum rule technique is employed to investigate pentaquark states with strangeness 
S = +1 and IJ n = 0i ,li ,0| ,l| . Throughout the calculation, emphasis is laid on the 
establishment of a valid Borel window, which corresponds to a region of the Borel mass, where the 
operator product expansion converges and the presumed ground state pole dominates the sum rules. 
Such a Borel window is achieved by constructing the sum rules from the difference of two independent 
correlators and by calculating the operator product expansion up to dimension 14. Furthermore, 
we discuss the possibility of the contamination of the sum rules by possible KN scattering states. 
As a result, we conclude that the 0| + state seems to be the most probable candidate for the 
experimentally observed + (154O), while we also obtain states with 0| , 
higher mass regions. 

PACS numbers: 12.38.Lg, 14.20.Pt 
Keywords: Pentaquark baryons, QCD sum rules 
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I. INTRODUCTION 

O + (1540) with strangeness S = +1 and baryon num- 
ber B = +1 is evidently a flavor exotic state with mini- 
mal quark content uudds [l[ . The first announcement of 
its experimental detection was made in 2003 and it 
has since confronted the hadron physics community with 
interesting novel phenomena and unanticipated problems 
that have not been solved until the present day. 

Presently, the experimental situation of O + (1540) 
seems to be rather unclear. After the CLAS collaboration 
has published several papers on their pentaquark search 
with high statistics where no signal of + could be 

found, many people now seem to believe that the pen- 
taquark does not exist after all and that the whole story 
was just "a curious episode in the history of science" Qj- 
There are, however still experiments that claim to ob- 
serve a sie nal of ©+(1540) HQ and therefore this issue 
should not be considered to be completely settled yet. 
Additional experimental results, which either unambigu- 
ously confirm the existence of O + (1540) or otherwise can 
eliminate it completely, are eagerly waited for. 

Theoretically, one not yet well understood property of 
+ (154O) is its unnaturally narrow width, which was re- 
ported to be even less than 1 MeV @ and which is very 
difficult to explain from our experience with ordinary 
baryons. Because + (154O) lies about 100 MeV above 
the KN threshold, one would expect the width to be 
much larger than the experimentally measured value. Of 
course, there have been many attempts to explain this 
narrow width of 6 + (1540) |10rll7l ]. but none of these ap- 
proaches has completely succeeded yet. 
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Another problem is the correct assignment of quan- 
tum numbers such as spin and parity to the ©+(1540) 
state. There are many studies, in which states with vari- 
ous quantum numbers were investigated using QCD sum 
rules [H-ii] or lattice QCD [UHl, but no consistent 
understanding has yet emerged. 

The main subject of the present paper is to deter- 
mine the quantum numbers (isospin, spin, and parity) of 
O + (1540) from a QCD sum rule approach. We there- 
fore study and compare the sum rules of states with 
I J* =0^ , l| ,0| ,l| . From this comparison, we 
aim to determine which of the investigated quantum 
numbers is the one that most likely has to be assigned to 
the O + (1540) state. Furthermore, we also look for possi- 
ble excited states below 2 GeV, that may be found in fu- 
ture experiments. For these purposes we use an improved 
version of the QCD sum rule method, which has first 
been proposed in (27j . The basic idea of improvement is 
to use the difference of two correlators to construct the 
sum rule. The continuum part of the spectral function is 
significantly suppressed by this procedure, which there- 
fore helps to find a valid Borel window, whose existence is 
a necessary condition for obtaining reliable results within 
the QCD sum rule technique. Moreover, we calculate the 
operator product expansion (OPE) up to dimension 14, 
which is indispensable for a sufficient convergence of the 
expansion. 

The paper is organized as follows. In Sec. II, the 
formalism of QCD sum rules is briefly reviewed. The 
details of our method, including the interpolating fields 
employed and the implications of the improvement men- 
tioned in the last paragraph are then explained in Sec. 
III. In Sec. IV, the results of the analysis for the various 
quantum numbers are given in detail. These results are 
then discussed in Sec. V and the conclusion is given in 
Sec. VI. Finally, the appendix is devoted to the numerical 
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results of the OPE and to the details of the establishment 
of the Borcl window for the various sum rules. 



In this study, we employ the usual "pole + continuum" 
parametrization for the spectral function, which appears 
in the imaginary part of the correlator 



II. FORMALISM 



A. QCD sum rules 

In the QCD sum rule method [111 [3(| , we compute the 
two-point function of various operators. It is defined as 



11(g) = i J d 4 xe lqx (0\T[r](x)r}{0)}\0) 

-n^ 2 ) ^ + n 2 ( g 2 ), 



(i) 



where rj(x) is a spin ^ operator. IIi(q 2 ) is called the 
chiral-even, and n^q 2 ) the chiral-odd part. 

Furthermore, the two-point function of a spin | Rarita- 
Schwinger type operator rj^ (x) is defined as 



(2) 



U(q) = i J d 4 xe^(0|T[^(x)7^(0)]|0) 

= -S^PIxOt 2 ) ^ + n 2 ( g 2 )] + ... 



where the dots stand for other Lorentz structures than 
g^ u . In the present study, we will only need the terms 
containing g^. We use the same notation as for the spin 
\ case and denote IIi(g 2 ) as the chiral-even, and II2 (g 2 ) 
as the chiral-odd part. 

In the QCD sum rule approach, we use the analytic 
properties of these two-point functions to extract infor- 
mation of the physical states that couple to the operators 
T] or 77^. Concretely, the analyticity of Eqs.([T|) and al- 
lows one to write down the dispersion relation 

, 1 [°° , ImlLfs) , , , , , 

"■i(Q J = — / as 5 — h (subtraction terms), (3) 

7T Jo s-q 

for i = 1,2. This equation has the same form for both 
the spin \ and spin | case. To handle a possible di- 
vergence in the integral of the right-hand side, usually 
the subtracted dispersion relation is used. This means 
that subtraction terms have to be added to this equa- 
tion, which contribute mainly to the high-energy part of 
the spectral function. The significance of these terms 
for the low-energy region is hence small. Moreover, the 
subtraction terms will disappear when the Borel trans- 
formation is applied as they are polynomials of q 2 . We 
will thus omit them in the following. 

The imaginary part of the two-point function Imllj(s), 
which corresponds to the spectral function of rj or 77^, 
satisfies the following spectral conditions: 

lmili(s)>0, v^Imlli(s) - India (s) > 0. (4) 

Thus, the positivity condition only holds for the chiral- 
even part, while the spectral function obtained from the 
chiral-odd part is allowed to have negative values. 



Imn,( S ) = 7 r(A t ) 2 (5( S -m 2 5+ )+0( S - St/i )Imn^( S ). (5) 

The 8- function for the ground state pole is justified by 
the experimental results, which show that the width of 
O + (1540) is very narrow. The potential contribution of 
the KN scattering states, which are not included in the 
expression above will be discussed in Sec. MIDI 

Furthermore, to obtain consistent results, it is im- 
portant to choose an appropriate operator in the two- 
point function, whose spectral function resembles that 
of Eq.© as much as possible. In other words, the cho- 
sen operator should couple strongly to the ground state 
pole (if it exists), leading to a large value of the residue 
(Aj) and at the same time should only have a small over- 
lap with the continuum states, which are included in the 
parametrization of Eq.© only above the threshold pa- 
rameter 8th- In the present study, we will try to construct 
such an operator by considering linear combinations of 
two independent operators and then fixing the mixing 
angles so that the results are consistent with the ansatz 
of Eq.©. 

While the low-energy part of the spectral function be- 
low the threshold parameter sth is phenomenologically 
parametrized as in Eq.©, the left-hand side of Eq.© 
and the second term of Eq. <j5j) are calculated analytically 
using the OPE. The results of this calculation can be 
generally expressed as follows: 

5 00 « 

U? PE (q 2 ) = £ C 2j ( q r- j log(V) + £ , 
j=o i=i [q 1 

n 2 OPi V) = E c 2j+1 tf) 5 - j MV) + E 



Cll+2j 

AW" 



(6) 



Here, C, contain various quark and gluon condensates 
and numerical factors. 

The next step in the calculation is to apply the Borel 
transformation, which is defined as 



£a/[IT(<7 2 )] 



lim 



i-q 2 ) 



2\n+l 



5? n ^ 2) < (7) 



where M is the so-called Borel mass. There are several 
reasons for using this transformation: firstly, the high- 
energy continuum part of the spectral function and the 
higher-order terms in the OPE are suppressed by the 
factors e~ s l M and 1/n!, respectively. This considerably 
improves the accuracy of the sum rules. Secondly, as 
already mentioned above, the Borel transformation re- 
moves the subtraction terms in Eq. © and therefore elim- 
inates possible ambiguities originating from these terms. 
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Substituting Eq.© into the dispersion relation of 
Eq.©, and applying the Borel transformation, the fol- 
lowing expressions can be obtained: 



(Ai) 2 e- ro e+/ M2 

rsth „ 5 

=-/ dS e~^Y, c ^- j +Y: 

J ° 1=0 j=l 



(-iyc 10+2j 

r(j)(M2y-i 



fi(M,s th ), 



(8) 



rath , 5 °° 

/ d S e-*/ M2y £C23 + is 5 - 3 +J2 

J° 3=0 3=1 



(-l)^ 11+ 23 

r(j)(M 2 y-i 



f2(M,S th ). 



(9) 



From these equations, the expressions for m e + and (Ai) 2 
can be extracted straightforwardly: 



m 2 e+ (M, 8 t h) 



1 



9/i(M, 8 th ) 



fi(M,sth) d(-l/M 2 



(A 4 ) 2 -/ l (A/. St , l ) e m e + W s -)/^ 2 



(10) 



:ll) 



Notice that rriQ+ can be calculated independently either 
from the chiral-even term fi(M,sth) or from the chiral- 
odd term f%{M, Sth). In this study we will mainly use 
/i(M, sth) to calculate m e + and refer to f 2 {M, sth) only 
for determining the parity of the investigated state. 

Here, the Borel mass M and the threshold parameter 
Sth arc variable parameters, which allow us to obtain 
distinct sum rules in Eqs. (fl"U]) and (TTTj) for each chosen 
value of M and s t h- Comparing these different sum rules, 
it is possible to extract information on the shape of the 
investigated spectral function and on the physical states 
that contribute to the sum rules most strongly. We will 
discuss this issue in detail in the later sections. 



B. Parity projection 

The parity of the presumed ground state pole can not 
be determined from the sum rule of the chiral-even (or 
chiral-odd) part alone, as rj(x) or r]^{x) couple to both 
states with positive and negative parity regardless of their 
own intrinsic parity. To this end, we use the parity- 
projected sum rules [37| to obtain information on the 
parity of the investigated state. In this method, instead 
of Eqs.(fT]) or ©, the "old fashioned" Green function is 
considered in the rest frame (q — 0): 



n o/ (g ) 



d 4 xe tqx (0\e(x )r)(x)r](0)\0) 



q=0 



n° / (g ) 7 u + n^( (Zo ) 



°fl 



(12) 



Here, only the spin ^ case is shown. This expression leads 
then to two independent sum rules for states coupling 
to rj(x) with positive and negative parity respectively. 
These are given as 



(mi 



|A±|V 

t 

1 rio 



7m 2 



dq [Irrtf]^ (qo) ± ImH^ (go 



(13) 



where the intrinsic parity of the operator has been as- 
sumed to be positive. In the opposite case the signs of 
the right-hand side have to be switched, q^ is the thresh- 

1 /2 

old parameter corresponding to s th in Eq.©. 

We now have three sum rules, Eqs.(|8|),(|9|) and (fl3|) . 
which must in principle give the same results for mg+ . 
However, the OPE of the chiral-odd part Ii 2 has turned 
out to contain ambiguous terms in the first power of the 
strange quark mass m s , related to an infrared divergence 
originating in the perturbative treatment of m s . To cir- 
cumvent this problem, we will use only the sum rule of the 
chiral-even part of Eq.© (where the divergencies do not 
occur) to calculate the mass of the ground state. Mean- 
while, Eq.(fT3|) will be applied in the chiral limit (-j|^y = 
1 , m s = 0) in order to determine the parity of the ground 
state. 



C. Borel window 

It is well known since the QCD sum rule method has 
been formulated, that the condition of an existing Borel 
window provides an essential check of the accuracy of 
the method. We define the Borel window as the region 
of the Borel mass where the following two conditions are 
satisfied. 



t \r\OPE 

[ highest order terms 



< o.i. 



dse"l^Imn OP£; (s) 
> 0.5. 



/ 

Jo 



dse~l^Imn OPE (s) 



(14) 



(15) 



Eq. (fT4")) is a necessary condition for the OPE to converge. 
It gives a lower limit for the Borel mass because the 
higher-order terms are suppressed for larger M as can 
be seen in Eqs.dU) and ((!]). On the other hand, Eq.(fT5|) 
ensures that the low-energy part of the spectral function 
dominates the sum rules, and that contributions from 
high-energy states do not deteriorate the result. It is 
therefore necessary that the high-energy part of the spec- 
tral function is sufficiently suppressed. As the suppres- 
sion is stronger for smaller M, this condition gives an 
upper bound for the Borel mass. 

One may wonder what the rationale for the numbers 
on the right-hand side of Eqs. (jT4")) and (Tl~5|) is. These 
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numbers are in fact chosen quite reasonably, which can 
be understood from the following considerations. The 
main uncertainties in a QCD sum rule calculation in most 
cases originate from ambiguities of the vacuum conden- 
sates which often have error bars considerably larger than 
10 %. This justifies the usage of Eq. (fl"4| as choosing any 
much smaller number than 0.1 on the right-hand side of 
this condition would be meaningless. Moreover, in order 
for the low-energy states below s t h to contribute most 
strongly to the sum rules and that thus an inappropriate 
parametrization of the high-energy states do not intro- 
duce too large errors, the right-hand side of ([15)) is also 
a natural choice, as is known from experience with sum 
rules of other baryons and mesons J36|. Therefore, one 
can have some confidence that the errors coming from the 
neglected higher-order terms of the OPE and the possi- 
ble inaccurate description of the spectral function above 
Sth are under control if both conditions of Eqs. (fT4")) and 
(fT5j) are satisfied. 

Let us now discuss the difficulties of establishing a 
Borel window in pentaquark studies. As has been dis- 
cussed above and also in [25|, H3, [H[ , the existence of a 
valid Borel window is essential for obtaining reliable re- 
sults within the QCD sum rule technique. Nevertheless 
almost all earlier studies investigating pentaquark states 
with QCD sum rules fliU^ I2H ] did not consider this 
problem and therefore these results should not be seen 
as to be conclusive. 

The reason why all these studies have ignored this is- 
sue, is that, in fact, it is very difficult (if not impossi- 
ble) to establish a valid Borel window in the conventional 
QCD sum rules described so far in this paper. There are 
basically two difficulties. Firstly, the convergence of the 
OPE of the correlator of a five-quark operator is slower 
compared to the case of nonexotic baryons containing 
only three quarks. This can be understood from a simple 
argument: graphs containing quark loops are multiplied 
by a factor of (<^:) 2 ' 1 (n: number of quark loops), which 
comes from the integration of the momenta in each loop, 
and the graphs of higher-order terms, where some of the 
quark loops are cut, are thus enhanced [Hj]. Therefore 
one can expect that the OPE starts to converge only af- 
ter all the quark loops are being cut. For the pentaquark 
calculation, this happens only for terms of dimension 12 
(or higher) and one therefore needs to calculate the OPE 
at least up to dimension 12 to make sure that all the 
terms with a possible large contribution are included. 
This problem can in principle be solved if one calculates 
the OPE up to high enough orders (which is a tedious, 
but straightforward task, if one uses the vacuum satura- 
tion approximation) . 

The second difficulty is more severe: because of the 
high dimensionality of the interpolating field of the pen- 
taquark (15/2 compared to 9/2 for an ordinary three- 
quark baryon) , the high-energy part of the spectral func- 
tion well above the presumed + resonance is enhanced 
and therefore in many cases dominates the behavior of 
the whole spectral function. This makes it very difficult 



TABLE I: The quantum numbers of the ud diquarks used in 
this study. a,b, . . . are color indices and C = i^j stands for 
the charge conjugation matrix. 



Diquarks 


I 


J 


7T 










+ 


£abc(v-aCd b ) 










tabciuICj^-yadb) 





1 




tabduIC-y^db) 


1 


1 


+ 



to obtain a large enough pole contribution in Eq. (|T5|) for 
establishing a valid Borel window. It has thus been a very 
hard task to make a reliable prediction on the resonance 
6+(1540). 

As discussed in the next section, this problem can be 
solved by a modification of the standard QCD sum rules 
technique, which consists of using, instead of a single 
correlator, the difference of two independent correlators 
to construct the sum rules. This will be our strategy in 
this paper. 



III. DETAILS OF THE METHOD 
A. Interpolating fields 

To carry out QCD sum rule calculations, one first has 
to construct appropriate operators, which carry the de- 
sired quantum numbers. These operators should be cho- 
sen to couple strongly to the state that one wants to 
investigate, although this is not always a trivial task. 
The interpolating fields that we use in the present study 
are described in this section. Our general strategy is to 
assemble two independent operators for each quantum 
number and set up general interpolating fields by consid- 
ering linear combinations of them. 

All the operators used in this study are built from two 
ud diquarks and an s antiquark, so that their KN com- 
ponent on the operator level is as small as possible. We 
therefore hope that these operators only have a small 
overlap with the KN scattering states while they should 
couple strongly to the possible pentaquark resonance. 
For orientation, the properties of the employed ud di- 
quarks are given in Table |U 

1. The IJ P = Oi* and 1^ states 

For the isosinglet case with spin \ , we use the following 
two operators. The same ones were used in (27|, where 

the IJ P = 0^ states were investigated with a similar 
strategy as in this paper. 

r)l=°{x) =e cfg [e abc ul(x)C^ 5 d b (x)} 

T -T 

x [e de fu d {x)C-i^ b d e {x)}^-i b Cs g (x), 
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rj 2 (x) =e cfg [e abc u a (x)Cd b (x)} 

x [e def uJ{x)C-/^ 5 d e (x)]^Cs^(x). 



(17) 



Here, a, b, . . . are color indices, C is the charge conjuga- 
tion matrix and T stands for the transposition operation. 
These fields both have positive intrinsic parity. They are 
constructed from a scalar diquark, a vector diquark and 
an antistrange quark operator in the case of r\\ and from 
a pseudoscalar diquark, a vector diquark, and an antis- 
trange quark operator in the case of r\2- To project out 
the spin ^ component, both operators have been mul- 
tiplied by 7^ and i]\ is furthermore multiplied by 75 to 
get the correct intrinsic parity. It must be remembered 
that even though the intrinsic parity of these operators 
is positive, it can couple to both states with positive and 
negative parity. 

By introducing a mixing angle 0^ 2 , a general operator 

can be constructed from rj^ =0 (x) and rj 2 =0 (x): 



COS^/^V) 



sin^/^- ^). 



(18) 



This is the operator that will be used in the actual cal- 
culation. Here, we are in principle allowed to choose any 
value for the mixing angle. The strategy for determining 
this free parameter will be given in the following subsec- 
tion. 

For the isotriplet, we employ the following two opera- 
tors, which have a form similar to the isosinglet case: 



r/ 1 (x) =e c fg[e a bcU a (x)C"f 5 d b {x)} 

x [edefuJix^f.deix^Cs^x), 



(19) 



contributions from states with spin i and |. It can be 
shown that the terms proportional to g^ v receive only 
contributions from the spin | states. Therefore, if one 



considers only the two terms 

n^( 9 ) = -.g^[n 1 ( (Z 2 ) 4 + n 2 {q 2 )] 



(22) 



the spin ^ contributions will automatically be eliminated. 
Note that there is a minus sign on the right-hand side of 
Eg. (f2"2")l , which is a consequence of the properties of the 
Rarita-Schwingcr field. 

To study the isosinglet states we employ the following 
interpolating fields 



(?) =e c / s [e a b c Ua(x)Cj 5 d b (x)] 

x [^de!U^{x)C-i^ b d e {x)\Ci^ (x), 



vL"(x) =e cfg [e abc ul(x)Cd b (x)} 

x [^defUd( x ) c lfi75d e (x)]j 5 Cs^(x) 



(23) 



(24) 



These operators have the same structure as the ones with 
spin I [Eqs. p^)) and (fTT)) ]. The only difference is that 
7^75^ front of Cs T has been omitted here, which al- 
lows the operators to couple to spin | states and lets 
the intrinsic parity become positive. As above, a general 
operator is then constructed by a linear combination of 
ri{=° and r?|=°: 



^° 3/2 vl=°(x) 



(x). (25) 



j] 2 (x) =e cfg [e abc u a (x)Cd b (x)] 

x [£de/uI(a;)C*7 Al 4(a;)]7 AI 75CsJ(a;) 



(20) 



The notation is the same as before. The difference to the 
fields with / = is that instead of a vector diquark, we 
here have an axial-vector diquark, which carries isospin 
1 = 1. Moreover, note that both the operators are mul- 
tiplied by (—75) to obtain positive intrinsic parity for 
T][ =1 (x) and ??2 =1 (' T )- Analogously to the isosinglet case, 
we construct a general operator by introducing a mixing 
angle d\j 2 : 

rf=\x) = co S 9 1 1/2 r 1 I r 1 (x)+sm9 1 1/2 r ] t 1 (x). (21) 



This is the same kind of operator that has been used 
in our previous work [28j . In this paper we will merely 
restate the results that have been obtained there in order 
to compare them with the results from the other quantum 
numbers. 

Finally, for the isotriplet case, we will use the operators 
given below: 



( x ) =e cfg [e abc u a (x)C~f 5 d b (x)} 

x i e defUd(x)C'y fJ ,d e (x)]>y 5 CsT(x), 



V2,v ( x ) =ecfg[e a bcU a (x)Cd b (x)} 

x [e de fUd(x)C'Y fJ ,d e (x)}CsT(x). 



(26) 



(27) 



2. The IJ P = Of* ana 



1| states 



The construction of the operators with spin | can be 
done in a similar fashion as for spin i. There are, how- 
ever, some additional steps arising from the properties of 
the spin | Rarita-Schwinger type fields. These are, for 
instance, discussed in [281 |39| and we do not repeat the 
details here. We only state the result of how the spin 
I components can be extracted. In the case of Rarita- 
Schwingcr fields, the two-point function of Eq.@ gen- 
erally contains various different tensor structures, with 



The structure of these operators is almost the same as 
the ones with quantum numbers IJ P = 1^ ■ Here 
again, compared with Eqs. (fT9]) and (|20[) the matrices 
7^75 have been omitted in order to construct Rarita- 
Schwinger fields which couple to spin | states and to 
adjust the intrinsic parity to be positive. As in all the 
cases above, a general operator r/j^ 1 is constructed from 

Wilt an( i rfTu-i which will then be used to formulate the 
sum rules 

< =1 (a0 = cos^t^Oe) + smel/rfsfix). (28) 
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B. Determination of the Borel mass, threshold 
parameter and mixing angle 

The Borel mass M appears in the formulation of 
QCD sum rules when the Borel transformation is ap- 
plied in Eq.([7]), the threshold parameter s t h in the 
"pole+continuum" ansatz of Eq. ([5]) , and the mixing an- 
gles 9 j in the general expressions for the interpolating 
fields in Eas.(fTg ]) ,(|3T ]) .(|25 j) and {25]). In this subsection, 
our strategy of determining these parameters will be ex- 
plained. 

Let us first discuss the question of how the Borel mass 
M has to be determined. As mentioned in the last sec- 
tion, it first has to be checked whether one can establish 
a valid Borel window from the sum rules. If not, the 
sum rules will not work and it will not be possible to 
obtain any reliable results from them. If one is able to 
find a valid Borel window, M has to be chosen within its 
boundaries. As will be discussed below, when the sum 
rules "work well" , the dependence of the results on M 
should be small and therefore it will not strongly depend 
on the exact position on M inside of the Borel window. 

Next, our strategy of determining the threshold param- 
eter Sth will be explained. Assuming that the low-energy 
part of the spectral function is dominated by a narrow 
resonance pole, the values of the resonance mass [given 
in Eq. ([fT)]) ] and the residue [Eq.([TT])] should not strongly 
depend on M and Sth- This is easily understood when 
one considers the (ideal) case, when the spectral function 
is given by a single 8- function below sth- Rewriting the 
right-hand side of Eq. (fTU]) . we obtain 

a(-iV) dse- s ' M \mn{ S ) 

f tH dse- s ' M2 lmR{s) 

f J° 2 (29) 

/ dse- s/M slmll(s) 

Jo 

r s th 

/ dse- s/M ImII(s) 
Jo 

If ImII(s) is a simple <5-function specified as ImII(s) = 
ir(\) 2 8{s - to 2 ) + 9(s - 8 t h)Tmrf(s), then Eq.{29]) gives 
to 2 and does not depend on M and Sth- On the other 
hand, if ImII(s) is described by some continuous positive 
curve, which corresponds to the scattering states, Eq . ([!?§)) 
should be a rising curve, because of the weight factor 
e -s/M- _ wn } cn suppresses the part of the integral with 
large s values when M is small. Furthermore, Eq.([29|) 
should have an increasing value when sth is raised, as 
higher values of s will be included in the integral. 

Following the arguments above, it can be understood 
that the threshold parameter sth has to be chosen so 
that the dependence of the calculated resonance mass 
and its residue is smallest, because this corresponds to 
the case of the largest contribution of a narrow ground 
state pole to the spectral function. On the other hand, 
if no such value for s t h can be found, we can assume 



that the spectral function is dominated by the scattering 
states. 

We therefore set up the following two conditions, by 
which we determine sth (called in the following the con- 
ditions of pole domination): 

1) A sufficiently wide Borel window exists. 

2) TO0+ (M, Sth) should only weakly depend on the Borel 

mass M and on the threshold parameter s t h- 

Condition 1) is essential to obtain reliable results with 
the QCD sum rule method, while 2) follows from the 
discussion above. The problem that arises here, is how 
we should quantitatively define the "weak dependence" 
of condition 2). In other words, how "weak" should the 
dependence on M and Sth be that one can be unambigu- 
ously sure not just to observe scattering states? This 
important problem will be discussed in the part of the 
result section, which deals with KN scattering states. 

Finally, the mixing angle 9j has to be fixed. To do 
this, we repeat the analysis outlined above for various 
values of 9j and at the end choose the one for which the 
conditions 1) and 2) are best satisfied. This concludes 
our discussion about the determination of the different 
parameters that appear in the sum rules. 



C. Establishment of a valid Borel window 

As was pointed out in Sec. Ill C[ it has so far been 
very difficult to establish a valid Borel window in QCD 
sum rule studies of pentaquarks. We aim to solve this 
problem by a modification of the standard QCD sum 
rules technique [I?} . The idea is to use, instead of a single 
correlator, the difference of two independent correlators 
to construct the sum rules. By this trick, it is hoped that 
we will achieve a large cancellation of the high-energy 
part of the spectral function, due to the restored chiral 
symmetry in this region. We will then be able to obtain 
a large pole contribution, which (if the OPE is calculated 
up to a sufficiently high dimension) will make it possible 
to establish a valid Borel window. 

To illustrate this point more concretely, let us consider 
the difference of two independent correlators that have 
been constructed in the second part of this section. Wc 
take as an example the operators with quantum numbers 
IJ P = 0| .In the same way as it was shown in [28| for 
operators of spin |, the operators of Eqs. p^)) and (fTT)) 
can be expressed with the help of the operators £i and 
£2, which belong to different chiral multiplcts: 

a = - 2(u T R Cd R )[(u T L C^d R ) - (u T R C^d L )}^Cs T R 

- 2{u T L Cd L )[(u T L C lli d R ) - (ulC^d L )]^Csl, 

£ 2 =2{u T R Cd R )[{ulC lll d R ) ~ {ulC 7fi d L )]rC^l 

+ 2(u T L Cd L )[(ulC^d R ) - {ulC^dL^Csl. 

(30) 
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Here, the color indices are omitted for simplicity. £1 
belongs to the (3, 15) © (15, 3) multiplct of the chiral 
SU(3)z, <g> SU(3) fl group with 4(1) right-handed and 1(4) 
left-handed quarks, and £2 to the (8, 8) multiplet with 
3(2) right-handed and 2(3) left-handed quarks. 

Using these chiral operators, rj[ =0 and r]2 =0 are given 

as 



r/=° 



(31) 



Thus, the correlator of the general operator of Eq. ([T51) . 
denoted as II^ 2 (g 2 , 0$y 2 ) = (v I ~°V I ~ )i can be expressed 
as follows: 



n? /2 (<M 



/2) 



\ [(66) 
1 



(66 



jcos(2e° /2 )[(a6) + (66)] 

isin(2^ /2 )[(a6>-(66>]- 



(32) 



The first term of this expression does not depend on the 
mixing angle #$y 2 , but * s ex P ec ted to couple strongly to 
the high-energy continuum states because this term can 
have perturbative parts. On the other hand, compared 
with the first term, the coupling to the high-energy states 
of the other two terms is expected to be smaller, which 
can be understood from the following arguments. The 
perturbative part of the second of Eq. (|3"2")) term vanishes 
because £1 and £2 belong to different chiral multiplcts 
and therefore at least one nonperturbative quark con- 
densate related to chiral symmetry breaking is needed 
to connect them. As the perturbative term largely cou- 
ples to the high-energy states, their contributions will be 
suppressed in this term. Considering the third term, it 
is possible to cancel the leading perturbative terms with 
an appropriate normalization of £1 and £2- Note that 
we here have implicitly used the positivity condition of 
the spectral function. As is seen in Eq.Q, this assump- 
tion is not necessarily valid for the chiral-odd part, but 
we employ in this paper only the sum rule of the chiral- 
even part and the parity-projected sum rules, where the 
positivity condition holds. 

Therefore, by taking the difference of two correlators 
with different mixing angles 9^ 2 and 0$y 2 , ^ ne ^ rs * term 
in Eq. (f3"2")l will be eliminated, and a strong suppression of 
the high-energy continuum part can be obtained. It will 
thus become possible to establish a valid Borcl window. 
Writing this difference down, we get 



n? /2 (g 2 ,C 2 )-n? /2 (<7 



- sin(0° /2 



??; 2 ){c OS (0? /2 



?i/ 2 )K66 



- (66 



? 1/2 )[(66 



(66)]}, 



which will be used to formulate the sum rules. It is un- 
derstood that the factor sin(0!jy 2 — ^?/ 2 ) has no influence 
on the mass of the ground state, calculated in Eq. (fTU|) . 

'1/2 = f an d on ly keep 



We hence fix it at 0j , 2 — 



6*1/2 + 9\j2 = 4*1/2 as a ^ ree P aram eter, which will have to 
be determined by the conditions stated in the last sub- 
section. 



D. Possible contribution of KN scattering states 

We have in this paper several times mentioned the pos- 
sible influence of the KN scattering states to the sum 
rules. Generally, if such scattering states have the same 
quantum numbers as the interpolating fields, they may 
always contribute to the sum rules to a certain extent, 
so we have to find a way to distinguish them from nar- 
row pole states that we are really interested in. We have 
already mentioned in the discussion of Eq. (|29[) . that in 
the ideal case when only one narrow pole is present in 
the spectral function, the results of the sum rules should 
not depend on the Borel mass M and the threshold pa- 
rameter sth- When only scattering states contribute to 
the spectral function, this behavior should change. The 
nature and extent of this change will be illustrated in this 
section. 

Let us first consider how the contribution of the KN 
scattering states to the spectral function should look like. 
It is known that the KN interaction in the S = +1 chan- 
nel is weak and slightly repulsive for / = 0, while the re- 
pulsion is stronger for 1 = 1 [Z(| 5l[ ■ As an illustration, 
we will here use phase space as a first approximation of 
the KN spectral function, which thus corresponds more 
closely to the 1 = case. Nevertheless, the qualitative 
behavior of the results of this section does not strongly 
depend on the detailed form of the spectral function and 
can therefore be considered to be quite general. 

In the case of spin | states, the contribution of KN 
phase space to the spectral function can be expressed as 
follows: 



1 



p(s = q ) = -\Xkn\ Im<j (-*) j (27r)4 



M 



(34) 



75" 



(<? ~ P) 2 ~ m2 + * e p 2 — + ie 



-(-75) 



Here, Xkn is the coupling of the used interpolating field 
to the KN scattering states, m and M are the masses 
of the kaon and the nucleon, respectively. The 75's are 
necessary because of the p-wave coupling of our interpo- 
lating field to the nucleon and kaon fields. 

Going to the rest frame (q = 0), the integral can be 
easily calculated and one gets 



(33) 



Pill) 



1 

4^2 



\akn\ -. YlO&N 



M), (35) 
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where En is the energy of the nucleon, expressed as 



E 



N 



M 2 



m 



2<7o 



(36) 



Therefore, the spectral function for the chiral-even part 
has the form 



-imn 



KN 



Wo) - ^ll A ^l - 



M 2 



(37) 



which contains contributions from both positive and neg- 
ative parity states. 

For spin | states, similar considerations can be applied, 
although there are some complications coming from pro- 
jecting out the contributions of the spin J = | states 
from the correlator. In this case, the spectral function is 
expressed as 



p(.s = , 2 ) = i ^ k # 



lJ H) 

\ y 'J (2tt) 4 



{q-p)ii{<i-v)v fl+M 



{q-pf 



iep* 



M 2 



(38) 



and the J = | projection is most easily done by applying 
the projection operator 



3q 2 



-Jlp.lv 



1 



2 (iplv ~ Quip) 



(39) 

Then, only considering the terms proportional to and 
again going to the rest frame, one obtains 



Pill) = 



1 IV 12 ( E \ 
1r> 2 2 \ A KN\ 



M 2 yi< 



4q 



g»v{loE N + M) 



(40) 



from which finally the spectral function of the chiral-even 
part can be extracted: 



-imnf w (^) = 



10 2 2 l-^i^AT | 2 



E N {E 2 N -M 2 f/ 2 



(41) 

This expression again contains contributions from both 
positive and negative parity. 

Next, we compute the results that would be obtained 
by the QCD sum rules if only the KN scattering states 
contribute to the spectral function. This means that 
we calculate the quantity corresponding to Eq. fTUl) or 
Eq. ([29|) . where for II(s), we now use the expressions ob- 
tained above. The results for spin ^ and spin | are given 
in Fig. [U 

It is clearly seen that while the dependence on the 
Borel mass M is relatively weak, the results depend 
strongly on the threshold parameter s t h- This can in- 
tuitively be understood from the fact that the spectral 
function containing only the phase space contribution is 



a fastly growing function with increasing energy. There- 
fore, the high-energy regions below the threshold param- 
eter will dominate the integral of Eq. ([%!))) , which then 
leads to a behavior as seen in Fig. [TJ with a strong de- 
pendence on Sth- Furthermore, this dominance of the 
high-energy states will make it difficult to obtain a large 
value for the pole contribution and to establish a valid 
Borel window. 

Note that we have here assumed that Xkn and X' KN 
to be constants with dimension [GeV 5 ]. They could in 
principle also have a dependence on go such as const, x , 
which would result in Eqs. (|3"T|) and (|4"Tj) being multiplied 
by gQ°. We have checked this case, finding that the re- 
sults are altered quantitatively, namely that the mass is 
shifted upwards while the difference between the differ- 
ent threshold curves increases. Nevertheless, our quanti- 
tative discussion above does not have to be changed, as 
our observation of a large dependence of the results on 
s t h is even more apparent in this case. 

The results of this section show that the dependence of 
Wq+ {M, sth) on sth provides us with an indicator of how 
much the KN scattering states contribute to the sum 
rule: a linear dependence of the same (or larger) extent as 
in Fig. [1] suggests a strong contribution of the scattering 
states, while a significantly smaller dependence indicates 
that a narrow pole exists and is the dominant structure 
in the spectral function. Similar arguments have already 
been discussed earlier in 14a. 



IV. RESULTS 
A. Preliminaries 

We summarize in this subsection general issues com- 
mon to the sum rules of all the various quantum num- 
bers, and explain the parameters, conventions and basic 
approximations used in the calculation. 

One important feature of the results for all quantum 
numbers is, that the perturbative term Co vanishes when 
the difference of the two correlators is taken. This cor- 
responds to the suppression of the contribution of the 
high-energy states as was discussed in the last section. 
The detailed results of the OPE of both the chiral-even 
and chiral-odd part are given in Appendix IA"1 

We will for all quantum numbers first investigate the 
sum rule for the chiral-even part and after that consider 
the parity-projected sum rules, where both the results 
of the chiral-even and chiral-odd parts have to be used. 
However, as mentioned before, the results of the OPE cal- 
culation of the chiral-odd part have turned out to contain 
some ambiguous terms in the first power of the strange 
quark mass m s , related to an infrared divergence. It is 
important to note here that this kind of divergence is ar- 
tificially arising because we are expanding our results in 
m s and arc ignoring higher-order terms. It should thus in 
principle be possible to remove this divergence by taking 
the full dependence on m s into account without resort- 
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1 1.2 1.4 1.6 1.8 2 1 1.2 1.4 1.6 1.8 2 

M [GeV] M [GeV] 

FIG. 1: The "mass of the ground state pole" of spin i (left) and | (right), obtained when only the KN scattering states 
contribute to the spectral function. Eq.([37[) and Eq. (|41[l have been used as the expression of the spectral function of the chiral 
even-part. 



TABLE II: Values of all the parameters used throughout this 
paper, given at a scale of 1 GeV [3r3 . |43| . The parameter 
k describes the possible breaking of the vacuum saturation 
approximation and is explained in this subsection. 

(qq) -(0.23 ±0.02 GeV) 3 
^4 0.8 ±0.2 

{q9 { ° q Gq) 0.8 ±0.1 GeV 2 
^■Gs) 0.8 ±0.1 GeV 2 

(ss) 

(if G 2 } 0.012 ± 0.004 GeV 4 
m a 0.12 ±0.06 GeV 

K 1 ~ 2 



ing to any expansion, although this may be difficult to 
achieve in practice. In any case, to avoid this kind of 
artificial ambiguity, the results of the parity-projected 
sum rules will be given only in the chiral limit (m s = 0, 
(ss) = (qq)). Hence, our strategy will be to use the 
results of the chiral-even part to calculate the mass of 
the investigated state, while we will employ the parity- 
projected sum rules to determine only the parity of the 
state. 

The values of the mixing angles and threshold parame- 
ters are obtained using the conditions of pole domination 
of Sec. MI B[ We will use the same values for both the 
sum rule of the chiral-even part and the parity-projected 
sum rules. 

The values of the condensates and other used parame- 
ters are given in Table [TTJ These are standard values for 
QCD sum rule calculations [36|,[43|, but they of course all 
have a certain range and the results will therefore depend 
on what exact values have been chosen for the conden- 
sates and other parameters. In the last part of this result 
section we will show to what extent the results will be 
influenced by the uncertainties of these parameters. 

Finally, re, the last parameter of Table HH will now be 



explained. It parametrizes the possible violation of the 
vacuum saturation approximation and is used as follows: 

(qqqq) = i^(qq) 2 , 
(qqqqqq) = K 2 (qq) 3 , 

(qqqga ■ Gq) = K(qq)(qg<j ■ Gq), (42) 
(qga ■ Gqqga ■ Gq) = K(qga ■ Gq) 2 , 

All the results shown below are obtained with k = 1, 
which means that the vacuum saturation approximation 
has been assumed. This approximation has been shown 
to be valid in the leading order of the large N c expansion, 
even though the w- corrections may be quite large. Wc 
have checked to what extent the results change when this 
approximation is broken up to values of k = 2. These 
changes will be included in the estimation of the error. 

The results of the various sum rules are given in the 
following. To allow a direct comparison between the dif- 
ferent quantum numbers, all the plots corresponding to 
the same quantity are shown in the same figure. 



B. IJ w = 0± ± 

We first show our obtained results for the isosin- 
glet, spin i case. This quantum number has been al- 
ready frequently investigated as a pos sible assignment 
for e+(1540) in QCD sum rules [Sill H, 13 and lat- 
tice QCD [29h32| . As for QCD sum rules, most of these 
calculations have problems in the establishment of the 
Borel window. In [271 ] , this problem was avoided by tak- 
ing the similar approach with that of the present paper, 
but we found some mistakes in the computation of the 
OPE beyond dimension 8. After correcting them, we got 
the present result which excludes the positive parity state 
obtained in (27| . 
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1 1.2 1.4 1.6 1.8 2 1 1.2 1.4 1.6 1.8 2 

M [GeV] M [GeV] 

FIG. 2: The mass of the pentaquark with quantum numbers 1,1 = 0|, l|, 0|, l| as a function of the Borel mass M, obtained 
from the sum rule of the chiral-even part. The arrows indicate the boundary of the Borel window for the middle value of the 
threshold parameter s t h- 



1. Sum rule for the chiral-even part 

Using the operators of Eqs.(|16 p .(|17 p . we first determine 
the mixing angle 4>y 2 ana - the threshold parameter s t h- 
The values that we have obtained are, <f>y 2 = —0.22 and 
y/sth = 2.2 GeV. Furthermore, checking the convergence 
of the OPE and investigating the value of the pole con- 
tribution, we have confirmed that a Borel window exists 
for 1.2 GeV < M < 1.6 GeV (for details, consult Figs. [3 
1 and 1 of Appendix [Bj . 

The calculated value of the ground state mass 
m@+(M, s t h) of Eq. (fl0]) is given in Fig. [2] (top left) as 
a function of the Borel mass M. The boundary of the 
Borel window for the case of y/sth = 2.2 GeV are indi- 
cated by the two arrows. One can see that the obtained 
value is about 1.5 GeV within the Borel window. Even 
though we have found a wide Borel window, the curves 
shown in Fig. [5] exhibit quite a large dependence on M 
and Sth, which suggests that the spectral function only 
contains KN scattering states and not a narrow pole. 
On the other hand, as will be shown later, the result of 
the parity-projected sum rules are fairly stable against 
M and Sth, which rather points to a narrow pole in the 



ground state. The interpretation these different results 
will be discussed below. 

The result of the residue, calculated from Eq. (fTTj) is 
given in Fig. [3J As well as for the mass, the results for 
the residue depend on M and Sth quite strongly. 



2. Parity-projected sum rules 

As already mentioned before, we will use the parity- 
projected sum rules in the chiral limit. This is justified, 
as we have confirmed in the sum rules of the chiral-even 
part that the qualitative behavior of the results does not 
change when this limit is taken. To show the strength of 
the contribution of the positive and negative parity states 
in the spectral function of the sum rule, the residues 
of the parity-projected sum rules [Eq. (|13|) ] are given in 
Fig. [2J It is clear from this figure that the negative 
parity states dominate and that therefore negative par- 
ity has to be assigned the state investigated in the last 
section. Furthermore, the mass calculated from the neg- 
ative parity sum rule of Eq. (|13|) is shown in Fig. [5J to- 
gether with the Borel window for the threshold parameter 
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FIG. 3: The residue (Aj ;1 ) 2 for the quantum numbers I J = Oh, l|, 0|, l|, obtained from Eq.JTTJ) for the sum of the chiral-even 
part. The value is given for three different threshold parameters. 



^Jsth = 2.2 GeV. As is seen in the figure, a valid Borcl 
window is established around 1.2 GeV < M < 1.3 GeV 
and the obtained value is consistent with the one of the 
chiral-even sum rule. Moreover, the dependencies on 
both M and Sth are very small, which in contrast to the 
chiral-even case rather points to a narrow ground state 
pole and not to KN scattering states. 

It is puzzling why the behavior of these two sum rules 
is so different, even though the contribution of the posi- 
tive parity states is very small, as shown in Fig. [4] Nu- 
merically, this can be understood from the fact that the 
chiral-even part is multiplied by an additional power of 
qo in the parity-projected sum rules [compare Eqs.([T]) 
and (p~2|) ] . which considerably changes the behavior of 
the sum rules in this case. Moreover, we have confirmed 
that even though the residue for positive parity state is 
small, it numerically has a large influence on Eq.® for 
the low Borel mass region. To illustrate this point, the 
contribution of positive and and negative parts, calcu- 
lated from Eq. (fTB"| in the chiral limit, are shown in Fig. 
|B1 The negative parity part clearly shows an unphysical 
behavior as it is almost constant, while it should be an 
exponentially increasing function in the case of a narrow 
ground state pole dominating the sum rules. Neverthe- 
less, around 1.2 GeV, its contribution is comparable to 



positive parity part and therefore has a strong influence 
on the result of the chiral-even part. Thus, the most rea- 
sonable explanation for these different results seems to 
be that the positive parity KN scattering states are con- 
taminating the results of the chiral-even part and there- 
fore lead to a large dependence on M and especially on 
Sth- We hence conclude that we have found some real 
evidence for a narrow ground state pole with IJ* = 0^ 
even though the situation is more ambiguous than in the 
other channels. 



C. IJ" = if* 

Next, the isotriplet, spin \ states are studied. As no 
isospin partners of + (154O) have so far been found ex- 
perimentally, it is currently believed to be an isosinglet 
state, but this assignment is not conclusive yet. Further- 
more, even if O + (1540) is an isosinglet state, a different 
isotriplet pentaquark state could exist at higher energies. 
We thus consider this state in the following paragraphs. 
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residues of the positive and negative parity sum rules for the quantum numbers IJ = 0^, l|, 0|, l|, obtained 
The curves are calculated in the chiral limit (m s = 0, (ss) = (qq)). 



Sum rule for the chiral- even part 



2. Parity-projected sum rules 



The method is essentially parallel to the isosinglct case, 
the difference being only that we employ the operators 
of Eqs.CEU) and ([H instead of Eqs.([T6|) and CEZ]). The 
values of the mixing angle 4>\^ 2 an d the threshold pa- 
rameter sth have turned out to be 4>\/ 2 ~ ~ 0.079 and 
v^T = 2.5 GeV. 



The mass calculated from Ea. ([TU)) is shown in Fig. [2] 
(top right), as before with the Borel window for the mid- 
dle value of the threshold parameter, indicated by the 
two arrows. The obtained value is about 1.6 GeV. Com- 
pared to the isosinglet spin \ case of Fig. [31 it is obvious 
that the dependence on the Borel mass M and especially 
on the threshold parameter s t h is small, which is positive 
evidence for a narrow ground state pole in the spectral 
function. The residues for the three different threshold 
parameters are given in Fig. [31 where we again get only 
a similarly mild dependence on M and Sth- 



We will follow the same method as in the isosinglet case 
and calculate the parity-projected sum rules in the chiral 
limit. As is shown in Table IIIIl the result of the chiral- 
even sum rule did depend on the strange quark mass 
m s quite strongly and one thus may wonder whether the 
procedure of taking the chiral limit is justified. But, as 
we will use this sum rule only to determine the parity of 
the state, we think that it is accurate enough to provide 
reliable information, because even though the mass value 
of the state may quantitatively change, it is improbable 
that the parity of the state will switch when this limit is 
taken. 

The residues of the positive and negative parity sum 
rules are compared in Fig. 31 In this figure, it is seen 
that both residues are similar in magnitude. (Note that, 
as we have taken the difference of two correlators, the 
residue (A]y 2 -J 2 can become negative. States with nega- 
tive residues can thus not be ruled out as unphysical like 
in the ordinary QCD sum rules with just one correlator.) 
On the other hand, the residue of the positive parity state 
is very unstable against the variation of the Borel mass, 
which suggests that it does not correspond to a narrow 



13 



2 
1.9 
1.8 
1.7 
1.6 
1.5 
1.4 
1.3 
1.2 

1.8 
1.7 
1.6 
1.5 
1.4 
1.3 
1.2 
1.1 
1 



IT = 0,1/2" 



1.2 



IJ = 0,3/2 



: 2.4 GeV 
: 2.2 GeV 
: 2.0 GeV 



1.4 1.6 
M [GeV] 



1.: 



: 2.0 GeV 
: 1.8 GeV 



1.2 1.4 
M [GeV] 



1.6 



1.8 




> 

■u 

p 



> 

O 



2.1 
2 
1.9 
1.! 
1.7 
1.6 
1.5 
1.4 
1.3 

2 
1.9 
1.8 
1.7 
1.6 
1.5 
1.4 
1.3 
1.2 



• 1,1/2" 




1.2 



1.4 



IJ = 1,3/2 



: 2.7 GeV 
: 2.5 GeV 
: 2.3 GeV 



1.6 1.8 
M [GeV] 



2.2 



: 2.4 GeV 
: 2.2 GeV 
: 2.0 GeV 



1.2 



1.4 1.6 
M [GeV] 



1.: 



2.4 




o| + , if 



3 + 



as a function of the Borel mass M. 



FIG. 5: The mass of the pentaquark for IF = 0| , IF = l\ , IF - 
The arrows indicate the boundary of the Borel window for the middle value of the threshold parameter sth- The curves are 
calculated in the chiral limit. 
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FIG. 6: The contribution of positive and negative parts of 
Eq.® for IJ = 0i, obtained from Ea. lfTg)) in the chiral limit. 
The value of the threshold parameter is ^fsth = 2.2 GeV. 



As expected form the result of the residues, the calculated 
mass of positive parity strongly depends on the Borel 
mass M and no stable region is found. In contrast, the 
results for negative parity are stable and consistent with 
the value obtained from the chiral-even sum rule. We 
therefore conclude that the parity of the state is negative. 
The mass values for the negative parity case are shown 
in Fig. [5] 

Compared to all other cases studied in this paper, the 
Borel window here seems to be unnaturally large. The 
reason for this is that the same phenomenon as in the 
upper left part of Fig. |H] has occurred, meaning that 
due to some cancellation in the integral of the spectral 
function above sth, a peak has emerged in the function of 
the pole contribution, which shifts the upper boundary 
of the Borel window to a high value and therefore leads 
to this very large Borel window. 



ground state pole that we are looking for. Meanwhile, 
the residue with negative parity is fairly stable and thus 
seems to be consistent with the chiral-even sum rule. 

The masses of both parity states are also calculated. 



D. LF = 



This quantum number has been already investigated in 
detail by the present authors in a recent paper [28j . We 
will not repeat the analysis given there and only restate 
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the most important results. 

The same strategy as in this paper was followed, mean- 
ing that the difference of two correlators was taken, and 
the values of the mixing angle and threshold parame- 
ter were determined from the conditions of pole dom- 
inance. The obtained values are <j>®/ 2 = 0.063 and 
y/sth = 2.0 GeV. This then leads to the mass values 
shown in Fig. [5] (bottom left), calculated from the chiral- 
even sum rule. The obtained value lies at about 1.4 GeV. 
The result shows both a small dependence on M and Sth, 
which suggests that a narrow ground state pole exists in 
the spectral function of this quantum number. 

The parity of the state is determined with the parity- 
projected sum rules, leading to Figs. @] and \5\ Fig. 3] 
shows that the pole strength is dominated by the residue 
of the positive parity state. Fig. [5] then confirms that 
the positive parity sum rules give stable results, which 
are consistent with the ones obtained from the chiral- 
even sum rule. 



E. IJ W = if* 

The existence of states with quantum numbers IJ P = 
1| have been suggested for instance by studies using 
the quark model [44| and the chiral unitary approach 
[45j . We further investigate them here using the QCD 
sum rule method. 



1. Sum rule for the chiral-even part 

The operators used are given in Eqs. (f2l>| and ([27} and 
the following values have been obtained for the mixing 
angle and the threshold parameter: 4>\i 2 = 0.024 and 
y/s^= 2.2 GeV. 

The results for the mass are shown in Fig. [3] (bot- 
tom right) together with the Borel window for ,/sth = 
2.2 GeV. As can be read off from the figure, a value 
around 1.6 GeV is obtained for the mass of the state. 
The dependence of the result on both M and s t h is weak, 
which suggests that a narrow pole is present in the spec- 
trum. 

The value of the residue (Xy 2 i) 2 i s given in Fig. [31 
where again only a small dependence on M and s t h is 
observed. 



2. Parity-projected sum rules 

We have obtained a consistent result for the positive 
parity channel, while no state below 2.0 GeV was found 
with negative parity. The two residues are shown in Fig. 
[3J where one can see that the magnitude of the positive 
parity residue is larger that the one of negative parity 
and that it is an almost completely stable against the 



variation of M. The calculated mass of the positive par- 
ity sum rule, shown in Fig. [5J moreover gives similar 
values as obtained in the chiral-even case, which do not 
strongly depend on M and Sth- We therefore conclude 
that positive parity has to be assigned to the investigated 
state. 



F. Estimation of the theoretical ambiguity 

As the last point, we have to investigate the depen- 
dencies of the results on the various parameters of Table 
HIl in order to obtain a quantitative estimate of the error 
inherent in our results. We will here use only the results 
of the chiral-even part for this estimation. 

The contributions to the errors for the different quan- 
tum numbers are given in Table IIII1 For instance, con- 
sidering the IJ P = 0^ case, we have already seen from 
Fig. [21 that the dependence of the mass value on M or 
sth leads to an uncertainty of about ±0.1 GeV. Among 
the other parameters, the result depends most strongly 
on (qq), which gives an uncertainty of about ±0.15 GeV. 
Similarly, raising the breaking parameter of the vacuum 
saturation approximation to n = 2 leads to an increase 
of the mass of about 0.15 GeV. Similar considerations 
lead to all the error contributions for the other quantum 
numbers given in Table IIIII 

Assuming that the various errors are uncorrelated, the 
final error estimations are then obtained by taking the 
root of the sum of all squared errors Srrii and rounding 
up: 



combined error ~ y^^(<5mi) 2 . (43) 

Note that this is merely a rough estimation, as there 
are additional errors coming from the truncation of the 
OPE and possible radiative corrections, that have been 
neglected in the current calculation. 

V. DISCUSSION 

The details of the results for the various quantum num- 
bers have been presented in the last section. Putting ev- 
erything together, these results can be summarized as in 
Table El 

A number of comments have to be made here. First 
of all, the statement "no state found below 2.0 GeV" in 
Table IIVI means that either no valid Borel window could 
be found or that the results of the sum rules did strongly 
depend on M and Sth and that therefore no evidence for 
a narrow ground state pole could be found. Concerning 
this point, in the case of IJ W = 0^ , the results of the 
chiral-even sum rule and the parity-projected sum rule 
are to a certain extent contradictory and we therefore 
have to put a question mark behind this conclusion. Fur- 
thermore, having found no narrow state in our sum rule 
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TABLE III: Contributions of the uncertainties of all parameters appearing in the calculation to the final error. Only values 
larger than ±0.05 GeV are explicitly given. These values have been obtained using the sum rule of the chiral-even part. 



u p = 9il Li! °il i£ 



M 


±0.10 GeV 


±0.05 GeV 


~ 


~ 


8th 


±0.10 GeV 


~ 


±0.05 GeV 


±0.05 GeV 


(qq) 


±0.15 GeV 


±0.20 GeV 


~ 


±0.10 GeV 


(ss) 


~ 


±0.05 GeV 


~ 


~ 


(qgcr-Gq) 

(qq) 


±0.05 GeV 


~ 


±0.10 GeV 


±0.05 GeV 


(sger-Gs) 
{ss) 


~ 


~ 


~ 


~ 




±0.05 GeV 


±0.20 GeV 


~ 


~ 


m a 


±0.05 GeV 


±0.20 GeV 


~ 


±0.05 GeV 


K 


±0.15 GeV 


±0.20 GeV 


±0.05 GeV 


±0.10 GeV 


combined error 


±0.3 GeV 


±0.4 GeV 


±0.2 GeV 


±0.3 GeV 



TABLE IV: Summarized results for all quantum numbers that have been investigated. The allowed KN decay channels of 
the respective quantum numbers are indicated in brackets. The mass values quoted here are obtained for the sum rule of the 
chiral-even part. 





Parity 


± 




J = h 


7 = 


no state found below 2.0 GeV 
(KN P-wavo) 


1.5 ±0.3 GeV (?) 

{KN S-wave) 


7 = 1 


no state found below 2.0 GeV 
(KN P-wavc) 


1.6 ± 0.4 GeV 

{KN S-wave) 


J = l 


7 = 


1.4 ±0.2 GeV 

(KN P-wavc) 


no state found below 2.0 GeV 
(KN D-wave) 


7 = 1 


1.6 ±0.3 GeV 
(KN P-wavc) 


no state found below 2.0 GeV 
{KN D-wave) 



calculation does not necessarily mean that such a state 
does not exist. It could happen that the spectral function 
is dominated by the KN scattering states and that the 
narrow states that we are looking for only couple weakly 
to the interpolating field that we have used. Nevertheless, 
we have constructed general operators from linear com- 
binations of two independent local operators, and have 
investigated all possible mixing angles and therefore the 
nonobservation of states shown in Table HVl has been con- 
firmed for quite a wide range of independent operators. 

Next, we discuss the physical implications of the ob- 
tained results. A question that comes to one's mind when 
looking at Table HVl is. why we can not observe J 71 " = -| + 

states while we are seeing the ones with J 7 ' = | . These 
states are in some models considered to be spin-orbit 
partners [l(|, so if these models are realistic and con- 
sistent with QCD, we should be able to observe both of 
these states. There are at least two possible explanations 
for our obtained results. One explanation could be that 
the states with J 7 * = i + in fact exist, but their coupling 
to the used operators are too small and/or the KN scat- 
tering contribution is too large, so that a narrow peak 



structure cannot be extracted. Another possible inter- 
pretation of the missing J 71 = ^ + states could be that, 
the spin-orbit partners of the spin | states are not the 
ones with spin i but with spin |. This would mean that 
9 + (1540) is indeed a very exotic state, as in this case the 
uudd quarks have to form a spin 2 state, which would 
then couple to the remaining s. This is of course only a 
very speculative picture, but it would be interesting to 
test it by calculating pentaquark states with spin |. 

Another important point, that needs to be discussed, 
is the interpretation of our results on the J T = i states. 
Such a state was also found in a lattice study (conducted 
only for the isosinglet state), where a resonance state 
was isolated from the KN scattering states [31| . Our 
results (especially in the isosinglet case) are somewhat 
ambiguous, and the errors are large, so it is difficult to 
draw any definite conclusions. In any case, whether such 
states turn out to be real pentaquark resonances or not, 
they most possibly do not correspond to the observed 
G + (1540) state, because J 7 * = \ states can decay into 
KN by an S-wave, for which the width is expected to be 
much larger than the observed value for + , which is less 
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than 1 MeV Q. Of course, in principle there may exist 
some so far unknown mechanism, which suppresses the 
width strongly and which would allow to assign the = 
2 quantum numbers to the + , but with our present 
knowledge and experience, this seems to be unlikely. 

VI. CONCLUSION 

We conclude from our results summarized in Table IIVI 
that the most probable quantum number candidate for 
6 + (1540) is IJ'* = 0| + . We have also found evidence 
for an isotriplet state 1| and two states with spin g 

(Otj and \\ ) at slightly higher energy. 

To obtain these results, we have employed the QCD 
sum rule method, whose reliability is improved by an- 
alyzing the difference of two independent correlators, 
by which the contribution of the high-energy continuum 
states is suppressed. Furthermore, by calculating the 
OPE up to dimension 14 it is made sure that the ex- 
pansion is converging well, and a valid Borcl window can 
be established. 

Considering the spin | states, although we could ob- 
serve some evidence for resonance states with IJ V = 0^ 

and lh in the region of 1.5 GeV, we have pointed out 
that the situation concerning the KN scattering states 
does not seem to be very clear and our predictive power is 
quantitatively very limited. Furthermore, as discussed in 
the previous section, we do not believe that these states 
correspond to O + (1540), because their width is expected 
to be too large to be consistent with the experimental 
value. 

Looking at the states with IJ n = 0^ and 1| , in 
both cases the values of the masses and residues show 
only a weak dependence on the Borel mass M and the 
threshold parameter sth- For the isosinglet case this was 
already pointed out in [2f|. This suggests that we are 
really observing narrow resonance states in the spectral 
functions of these quantum numbers. As no isospin part- 



ners of the O + (1540) have so far been found, it is believed 
to be an isosinglet, which leads to our conclusion that the 
O + (1540)is likely to be a state with quantum numbers 
IJ V = 0| + . The isotriplet state l| + is predicted to 
exist somewhere above the isosinglet, so it may be in- 
teresting for future experiments to look for this state. 
One nevertheless has to be cautious when interpreting 
the current results, as we cannot make any real quanti- 
tative prediction about the width of the state with the 
present method. Therefore it is difficult to say whether 
the predicted isotriplet state is narrow enough to be un- 
ambiguously detected in an experiment. 

The width of a state can be obtained from the QCD 
sum rule technique by calculating three-point functions 
of appropriate currents, and it would be interesting to see 
whether it is possible to obtain a value consistent with 
experiment for the 0| state and whether the IJ^ = 
l| state is really narrow enough to be experimentally 
observed. Furthermore, it is important to check whether 
our conjecture of the large widths of the J n = A states 
is really true or not. These issues are left for further 
studies. 
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Appendix A: Results of the operator product expansion 

We obtain the following result for the OPE in terms of the parameters Ci defined in Eq.([6]). Note, that we here 
give the values of Q after the difference of the two correlators is taken, and that we have used @j — 6j = \ and 
(f)j = 9 j + 6 j . After showing the results of the chiral-even part (up to terms proportional to m s ), the chiral-odd part 
is given in the chiral limit. 

The used abbreviations arc G 2 = G a llv G a > lv and a ■ G = cr^^G^, the A a being the Cell-Mann matrices, g is the 

a 2 

coupling constant of QCD, giving a s = |— . The values of the condensates and the strange quark mass are shown in 
Table HI 



1. u p = o^ 

a. Chiral-even part 



n n n l \ 9 V G ' 1 ) ,o n (<7<?) 2 • ,o m s (sgcr ■ Gs) 
G o =0, C 4 = - 2 u 3 . 57t6 cos0 1/2 , G 6 = - ^jsin^j 2i%6 C0S ^l/2' 

(gg) (ggg ■ Gg) . 

8 = 2 9 3 2 7T 4 ( 7C ° S( ^l/2 + 34sln< ^l/2). 

(qgcr-Gq) 2 . (<7i?) 2 (if G 2 ) „ . 

10 = 2 13 3 2 tt 4 ( 22coS( ^i/2 + 299sm0 1/2 ) — (6cos0 1/2 + 61sm<^ /2 ) 



13m s (^G 2 )(sga-Gs) m s (qq) 2 



.s.s 



^14^4 C0S *l/2 - 2 3 3 2^ 3111 



in 0i/2 1 



_ 2(gg) 4 . (gg)(g g q-G g )(^G 2 ) ^ . 

C12 =— sm</2 + 2 12 3 3 tt 2 (65 cos <^ /2 + 418 sin <^ /2 ) 

m s (qq) 2 (sga-Gs) , . 7m s (qq)(ss)(qgcj ■ Gg) . 
2W ( 3cos ^i/2 + 13sm0 1/2 ) + sm0 1/2> 

G14 = + ^ 2 9 3 3 tt 2 (17 cos + 8sm<^ /2 ) 

, 19m s (s S )(q gt j-Gq) 2 . 19m a (gg) 2 (ss) (^G 2 ) . 
+ sm ^1/2 + 2^4 Sln V 2 ■ 

b. Chiral-odd part 

n n r, (* s ) ^ 7(sga-Gs) . ^ _ 5(ss)(^G 2 ) . . 

161(^G 2 )(s gt x-G s ) . (gg) 2 (^s) . 

9 = 2i53V Sm ^ 2 ~ 2 2 3 2 ^ 2 sm ^/2' 

7(qq) 2 (sga-Gs) 7(qq)(ss)(qgcr ■ Gg) 
C " = ^63^ ( cos ^i/2 + 2sm0 1/2 ) + sln 0i/ 2 , 

25(g g ) 2 (ss)(^G 2 ) . (qq)(qg<T-Gq)$g<T-G S ) , AriA ^ 1Q _ Q . 
Ci3 = ^3! sm^ /2 2 l2 3% 2 ^ (424cos<^ 1/2 + 1979 sin 1/2 ) 

19(g S )(ggg-G g ) 2 . 
2W Sm ^ 2 - 



2. /J^lf* 



a. Chiral-even part 



n n n (lf G2 ) ,1 ^ (??) 2 • ,1 m s (sga-Gs) x 
G o =0, C 4 = - 2 u 3-57r6 cos0 1/2! C 6 --^2^4sm0 1/2 — cos0 1/2 , 

8 = 293V (3cos0 1/2 + 44sm0 1/2 ), 

C10 = - ^^/ (IScos^ + 361sin^ /2 ) + M^l^l(18co S ^ /2 - 29sin# /2 ) 

13m s (^G 2 )(sga-G S ) . m s (qq) 2 (ss) . x 

2W C ° S ^ ~ 2 2 3% 2 Sm 01 /2 . 

r 2(gg) 4 . . (gg)(g,ga-G g )(^G 2 ) ^ oK . 

G12 = 33— sm</> 1/2 2 12 3 3 tt 2 (135cos0 1/2 - 158sm0 1/2 ) 
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, m s (qq) 2 (sga-Gs) , , 19m s (qq)(ss)(qga ■ Gq) . x 
+ ^3^2 (9cos<^ 1/2 + 13sm0 1/2 ) + sm^ 1/2 , 

31(qq) 3 {qga ■ Gq) x m s (gg) ( g gcr ■ Gg) (sgcr ■ Gs) x x 

Ci4 = ^ sm0 1/2 + 2103^2 (69cos<^ /2 + 122sm0 1/2 ) 

, llm^Xggg-Gg) 2 . ! , 23m a ( gg ) 2 (ss)(^G 2 ) . 
+ 2^3^ Sm</>1 / 2 + 2734 sin ^1/2- 



b. Chiral-odd part 



(ss) , 7(sg<r-Gs) , 5(ss)(^G 2 ) , 

Ci=0, G 3 = - 2ll3 2 57r6 sm0 1/2 , G 5 = + 2l43Vi sin0 1/2! G 7 = + ^f-g sin0 1/2 , 



161(^G 2 )(s 5( r-G S ) . . ( gg ) 2 



C9 = "21532^4 sin0 1/2 - 2 . 327r2 sin^ 1/2 , (A4) 

(qq) 2 (sg<J ■ Gs) x . x 19( gg ) (ss) (g.gcr ■ Gg) . x 



= 2 6 3 2 tt 2 ( 3cos ^i/2 + 19sin0i /2 ) + — — 2 5 327r " 2 -sin^i/2 

29(gg) 2 (ss)(^G 2 ) . (qq)(qg<T-Gg)(sgvGs) ^ x . x 

G 13 = sin </> 1/2 2 12 3 3 tt 2 ^ ^ / 2 ^/ 2 ' 

ll(s S )(g g a-Gg) 2 



2 7 3 2 7T 2 



3in^ /2 . 



3. IJ P = 01* 



See 



4. JJ P = If* 



a. Chiral-even part 



nan (lT G2 ) ,1 n ^<?) 2 . ,1 . m s (sga ■ Gs) . 

' 4 = 2^33 -5^ cos ^/2' Ce = 263 2 57r 4 Sm03 / 2 + 2 14 3 ■ 5tt6 cos ^/2, 

g»=- ^ > 2 ^ g \ co S ^ /2+ 54 S in^ 2 ), 

Cl ° = ^W )2 ( 15cos ^3/2 + 677sin^ /2 ) - ^fj 2 } (2cos^ /2 + llsin$ /2 ) 

13m,(^G 2 )(s g a-G S ) m s (gg) 2 (ss) . x 

+ ^1633^4 cos^ 3/2 + sm^ 3/2 , (A5) 

r (OT) 4 • ,1 , 5( g g)( g ffa-G g )(^G 2 ) ^ . 

G12 =-p-sin<?!>3 /2 + 2 Ti 3 2 ^ 2 ( 3cos ^3/2 + 2Osin0 3/2 ) 

m s (qq) 2 (sg(T ■ Gs) x . x llm s (qq)(ss)(qg<r ■ Gq) . t 

2W ( C ° S 03 / 2 + 10 8111 <M ^2^ 8111 ^3/2 . 

97 (qq) 3 (qga ■ Gq) x m s (qq)(qgc7 ■ Gq)(sgcr ■ Gs) x x 

14 = 2*3* sm ^3/2 2 11 3 3 7r 3 (23cos</» 3/2 + 18sm0 3/2 ) 

llm s (ss)( gg <7-G g ) 2 x 25m s (gg) 2 (ss)(^G 2 ) 

2 9 3 4 tt 2 03 / 2 2^ 03/2 ' 



19 



b. Chiral-odd part 

Ci=0, C 3 = -21333^6 8111 ^3/2 . C 5 = + 2l3337r6 sm0 3/2 , C 7 = 2143 3^ 4 sm0 3/2 , 

5(^G 2 )(?gq.fa) . (gg) 2 (^) . , x .... 

9 = + 2^3%4 8111 ^3/2 - 2 3 327r 2 Sm ^3/2 . (A6) 

{qq) 2 (sga-Gs} x . x 19(qq)(ss)(qg<7 ■ Gq) . , 
C H = ^2^2 ( COS< ^3/2 ~ 15sm0 3/2 ) + 2^5-5 S1110 3/2 , 

13 = 2635 sin ^3/2 + 2 n 337r 2 (17cos0 3/2 - 216sm0 3/2 ) 

ll(ss)(g gf x-G g ) 2 . j 

^ — ; — ^ Sill 05 /Q . 

2 6 3% 2 V3/2 

Appendix B: Establishment of a valid Borel window 

In this Appendix, we explicitly show that a Borel window has been obtained for the sum rules of the chiral-even 
part for the various quantum numbers. 

First, the convergence of the OPE is checked. This is done by calculating the left-hand side of Eq. (fl4")) . The results 
are given in Fig. [7] Additionally, the right-hand side of Eq. ([SJ added order by order is shown in Fig. |8] to get a 
better idea of the behavior of the expansion. Subsequently, the pole contribution of Eq. (|15[) is investigated. The 
corresponding plots are given in Fig. [9] for the various quantum numbers. 
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FIG. 8: Contributions of different dimensions to the right-hand side of Eq.([8]) for I J = 0|, I J = 1^, I J 
in succession. The expansion in most cases starts to converge after terms up to dimension 10 are included. 
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FIG. 9: The pole contribution to the sum rule with I J = 0± I J = 1±, I J = Of, I J = if in comparison with the contribution 
of the high-energy continuum states. These plots, together with the ones in Fig. [7] finally lead to the Borel windows referred 
to in the result section. 



